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Abstract. We present a method designed to estimate the noise power spectrum in the time domain for CMB experiments. 
The noise power spectrum is extracted from the time ordered data avoiding the contamination coming from sky signal and 
accounting the pixellisation of the signal and the projection of the noise when making intermediate sky projections. This 
method is simple to implement and relies on Monte-Carlo simulations, it runs on a simple desk computer. We also propose a 
trick for filtering data before making coadded maps in order to avoid ringing due to the presence of signal in the timelines. 
These algorithms were succesfully tested on Archeops data. 
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1. Introduction 

Measuring the Cosmic Microwave Background anisotropics 
angular power spectrum has proved to be a powerful cosmo- 
logical tool giving direct information on both the cosmological 
parameters and the primordial Universe through the origin of 
the initial perturbations. 

The usual method used nowadays for measuring the CMB 
temperature fluctuations is to scan the sky with a detector (with 
a bolometer for instance) that measures the temperature varia- 
tions in a given beam. This provides time data streams along 
with pointing directions that are subsequently compressed into 
a sky map. In the general case, the noise in the time streams is 
not white (due to atmospheric contamination, spurious noises, 
electronic and detector temperature fluctuations). The sky is 
also in general not regularly scanned so that at the end, the 
noise in the maps is neither homogeneous nor white. 

This causes difficulties when trying to extract the signal 
power spectrum from the maps as it has to be separated from 
the noise. Fully optimal methods finding the maximum like- 
lihood solution for both the map and the CMB power spec- 
trum have been proposed such as MADCAP (Borrill 1999 1 
but they require the inversion of large covariance matrices 
on parrallel supercomputers and are hard to implement for 
the present generation of experiments that cover a large por- 
tion of the sky and have a large number of pixels (WMAP, 
Archeops and Planck in the near future). Alternative meth- 
ods that replace the full inversion by a Monte-Carlo simu- 
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lation of the noise angular power spectrum have been pro- 
posed JHivon et al. 2 001 Sz apudi et al. 200 Ij l. The results ob- 
tained with such methods are satisfactory and they are little 
time consuming compared to the maximum likelihood meth- 
ods. The analysis process is done in the following way in these 
Monte-Carlo methods : 

1 . noise power spectrum estimation in the data: it is very 
important not to be contaminated by signal at this stage. 
This is particularly difficult when dealing with large por- 
tions of the sky where galactic dust and clouds can be 
found. We propose in section |2l a method that provides an 
unbiased estimate of the noise Fourier power spectrum by 
correcting an initial guess via Monte-Carlo simulations. 

2. Fast map making: can be a simple coaddition of filtered 
timelines or a optimal iterative mapmaking that converges 
to the maximum likelihood solution for instance with a con- 
jugate gradient method. This map making has to be very 
fast as it is repeated for each Monte-Carlo realisation. We 
propose in section|3a method for making simply coadded 
maps of filtered data avoiding the ringing due to the effect 
of the filtering on bright sources such as the Galactic plane. 

3. Ce estimation on the map: This is done very easily 
with the Healpix package JGorski et al. 19981 1. The angu- 
lar power spectrum obtained here is called pseudo-C( as it 
differs from the true signal Cc due to various effects: pres- 
ence of noise, beam and time domain filtering effects on 
the signal, incomplete sky coverage that degrades the reso- 
lution in harmonic space. All these eff"ects can be accounted 
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for and coiTected with the MASTER JHivon et al. 200 1> or 
SpICE (jSzapudi et al. 2001 1 algorithms. 

In this article we propose methods for solving pointsn](in sec- 
tion|5J and|2(in section|5}- These methods have been sucess- 
fuUy used for the Archeops data analysis (Ben oit et al. 2003al 
IBenoit et al. 2003b We used the MASTER algorithm for go- 
ing from pseudo-Ct to true Q and for estimating the error bars. 

2. Noise spectrum estimation 

We present in this section the method we developped for deter- 
mining the noise Fourier power spectrum of a time data stream 
composed from CMB fluctuations. Galactic dust and coloured 
gaussian and stationary noise. The data is denoted as an ele- 
ments vector d (N is the number of time samples) and is com- 
posed of noise n and sky signal s: 

d - s + n (1) 

The noise power spectrum is given by P{v) - {h*h) where ~ is 
the Fourier transform and {) the ensemble average over many 
realisations. The point here is to avoid contamination from sig- 
nal (CMB and in a larger part Galactic dust) that would tend to 
overestimate the estimated noise power spectrum and therefore 
N( leading to a biased C( spectrum at the end. 

2.1. Raw estimation of tine noise spectrum 

The first raw estimation of the noise power can come from the 
data itself (that is n ^ d) as we are dealing with low signal to 
noise experiments (at least in the time domain as the noise is 
reduced by redundancy in the map). This is a correct approx- 
imation only in the case where the observed region of the sky 
is free from any bright source. In most recent experiments such 
as WMAP and Archeops (as well as in Planck) the instrument 
covers a large portion of the sky by making large circles that 
cross the Galactic plane twice per rotation. There is therefore 
a very bright source in all the portions of the time streams that 
prevent us from using the naive approximation n ^ d. 

A simple method to get rid of the signal contribution in the 
noise power spectrum estimation is to project the data onto a 
map of the sky portion observed. We apply a prior highpass 
to the data to select the frequancies we are interested in and to 
get a realistic map. In such a map, the noise level is reduced but 
the signal is almost unchanged. Reading back this map with the 
scanning strategy leads to a timeline where the signal to noise is 
greatly improved. This latter timeline can then be subtracted to 
the initial giving a noise estimate that is, in principle, free from 
signal contamination. In the matrix notation that has become 
common in CMB analysis this corresponds to the following 
operation: 

n,an.^d-A{A^Ay' A^d (2) 

where A is the pointing matrix. The noise spectrum estimation 
is then just obtained from the Fourier transform of riraw 



2.2. Biasing effects 

If one rewrites equation (|3 using the definition Q, one gets: 

n™„ = « - A{A^Ay^A'^n + (^s - A {a^ Af^ A'^ sj 

= n - r„ + rj 

where r„ and can be interpreted in the following way: 

residual noise (r„): The number of time samples per pixel is 
finite (in Archeops and for a resolution of 15 arcminutes we 
have in average 25 samples per pixel with a very asymetric 
distribution) and therefore there remains some noise in the 
projected map. This noise in the map, once re-read with the 
scanning stategy is strongly correlated to the initial noise in 
the timeline. We therefore slightly underestimate the noise 
spectrum. Finally r„ increases if the pixel size is reduced. 
This residual noise is also discussed in |Stompor et al. 2002| 

residual signal (r,): Because of the pixellisation of the pro- 
jected map, the high temporal frequencies of the signal 
(mainly dust a low galactic latitude or point sources) are 
badly subtracted from the initial time stream. This leads to 
an overestimation of the noise power spectrum. Finally 
decreases if the pixel size is reduced as it only originates 
from the large size of the pixels compared to the frequen- 
cies present in the data. 

When computing the noise estimate power spectrum, one gets: 

Prawiv) = {(n-fn+f,T(n-h+f.s)) (4) 

= {(h-r„y (n -?„)) + {f*f,) (5) 

because noise and signal are supposed uncorrected. The first 
term can be rewritten as (1 - e(v))P(v). Simulations show that 
e(v) does not depend significantly on the noise power spectrum 
in the frequency range considered here (v > 0.3Hz), but essen- 
tially on the scanning strategy. The second term is additive and 
depends on the signal contained in the map and can be defined 
as e.s(v). One can therefore rewrite equation (|5} as: 

P™.-(v) = (l-e(v))P(y) + e,(v) (6) 

In the case of the Archeops experiment these two effects are 
not negligible at the resolutions we are interested in. As they 
cannot be reduced simultaneously by increasing or decreasing 
the pixel size, both have to be accounted for in order to cor- 
rectly estimate the noise power spectrum. In the following, we 
have chosen to use relatively small pixel size (15 arcminutes) 
in order to reduce the contribution of to the bias as its de- 
termination, as will be seen later, relies on calibration ans is 
model-dependent. 

If one can determine the shape of both e(v) and ej(v), then 
debiasing the noise power is straighforward by inverting equa- 
tion (|6}. 

2.3. Estimating the noise effect 

The simplest way to estimate e(v) is to make noise only Monte- 
Carlo realizations of the data streams with spectrum Praw (this 
assumes that Prawiv) is close enough to P so that e is computed 
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Fig. 1. Noise bias e(v) (multiplicative). 



Fig. 3. Signal Bias es(v) (additive). 
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Fig. 2. Power spectrum of the noise bias e(v)xf( v) in harmonic Fig. 4. Power spectrum of the signal bias es{v) in harmonic 
space. space. 



with enough accuracy altough the algorithm could be iterated) 
and perform the operation of equation Q on each realisation. 
The average power spectrum of all realizations is biased by 1 - 
e(v) that can be calculated by taking the ratio of the input power 
spectrum to the output. The result obtained is shown in Fig.^ 
Large oscillations can be seen at lower frequencies with peaks 
at harmonics of the spinning frquency. This is not surprising as 
we expect the noise to project onto the sky at these frequencies. 
It is therefore at these particular frequencies that the noise bias 
is the largest, leading to large values of e(v). The noise bias 
transformed into harmonic power spectrum using the Archeops 
pointing strategy is shown in Fig. |2l showing that this effect 
should not be neglected as it amounts to a significant fraction 
of the error bars (shown in black in Fig.|6j. 

2.4. Estimating the signal effect 

The signal bias estimation does not require averaging over re- 
alizations as it depends only on the sky signal, dominated by 
Galactic dust. The estimation is simply done by adding real- 
istic Galactic dust signal on a noise only timeline (with spec- 
trum Praw(v)). The realistic simulated Galactic dust signal is 



obtained from |Schlegel et al. 1998| dust maps filtered with the 
corresponding time constants, beam and electronic filter. The 
noise-Hsignal timeline is then deglitched and filtered the same 
way as for the real data. We then remove the initial noise before 
applying the operation of equation (|2j. The galactic bias ej(v) 
is the Fourier power spectrum of the residual signal. The sig- 
nal bias €s(v) is shown in Fig.|3lin Fourier space and in Fig.|3 
Again, the power spectrum bias amounts to a significant frac- 
tion of the error bars, showing that the signal bias should not be 
neglected. The level of galactic bias highly depends on the level 
of Galactic signal in the data and therefore the calibration. The 
biggest uncertainty in the estimation of ^^(v) therefore arises 
from the uncertainty in the calibration typically around 10% 
(in temperature). It is finally relatively small compared to the 
signal as it only applies to the signal residual that is already a 
few percent of the whole signal. 

2.5. Accuracy on the reconstructed noise spectrum 

After debiasing of both effects, we obtain a noise spectrum that 
should, in principle reflect the initial noise spectrum very accu- 
rately. In order to estimate this effect, we performed 60 Simula- 
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Fig. 5. Results of the noise spectrum reconstruction obtained with realistic realizations. The top panel shows the initial noise 
spectrum in green and the average reconstruction over the 60 realisations in black. The lower panel shows the relative residual 
bias which is much smaller than 2 percent at all frequencies. 



tions of the debiasing procedure based on the same inital noise 
spectrum (a realistic Archeops noise spectrum). The result is 
shown in Fig. |5l The noise reconstruction is better than two 
percents at all frequencies which is sufficient for accurate noise 
spectrum subtraction with MASTER in the case studied here. 



3. Making maps with filtered data : filtering without 
ringing 



Making coadded maps using heavUy filtered data (in the 
Archeops case only frequencies between 0.3 and 45 Hz are 
kept) is a simple and accurate alternative to optimal mapmak- 
ing when the knee frequency of the noise is low enough to pre- 
vent the filtering from removing too much signal in the low 
frequency part. Even when this is the case, one has to make 
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sure that the filtering (obtained by setting low Fourier modes to 
zero) has the expected effect and does not pollute the data. 

Spurious effects from filtering do arise if there is strong 
signal in the data. A good example of such a strong signal is 
the Galactic plane crossings that occur twice per rotation with 
Archeops. As the galactic peak is essentially non stationary the 
filtering creates ringing signal around the galactic plane in the 
maps. The structures created are not only located just around 
the galactic plane (where they are the strongest) but all over 
the map at a level comparable to the CMB anisotropics that are 
searched for We therefore need a trick to remove these ringing 
effects, that is, filtering without ringing. 

The usual trick used to reduce ringing effects on filtered 
data is not to use too sharp filters. We used a filter with a sine 
shape going from to 1 in 0.06 Hz. We obtained this way good 
reduction of the ringing effect but we wanted to improve further 
the mapmaking. 

As the ringing arises from the presence of the Galaxy cross- 
ing in the data, removing the Galactic signal should reduce 
drastically the effect. A first step is to remove simulated galac- 
tic signal from |Schlegel et al. 1998| but the agreement between 
the simulated maps and the true Galactic signal is not perfect 
and there remains high peaks in the data. This method also in- 
troduces external data into the experimental data, a procedure 
which has to be avoided if possible. We therefore decided to 
remove an estimate of the Galactic signal done with the exper- 
imental data itself. This is done in the following way: 

- We first fit the data with a set of low frequency 
functions (sine functions multiplied by a gaus- 
sian ( |Bourrachot et al., in prep.^ ) not taking into account 
data located in bright regions of the Galaxy (using a mask 
made from Schleg el et al. 1998| galactic dust maps). The 
low frequency template obtained is then free from Galactic 
contamination. 

- We then remove this low frequency template from the ini- 
tial data and obtain a timeline containing mostly high fre- 
quency noise and Galactic signal and no low frequency 
stripes. 

- This timeline is then simply coadded to obtain a map. We 
then set to zero the pixels out of the Galactic mask used be- 
fore so that the map contains Galactic signal in the region 
known to be contaminated and zero elsewhere (in partic- 
ular, in the region where the CMB analysis is to be per- 
formed). 

- This map is read with the scanning strategy and the result- 
ing timeline is subtracted from the initial data. The resulting 
timeline is stricly identical to the initial one in the CMB re- 
gion and differs in the Galactic region just by the estimate 
of the Galactic signal. 

- This last timeline is then filtered using a sine shape filter 
between 0.3 and 45 Hz and coadded to make the map that 
will be used for the CMB analysis. 

This method allows a very efficient removal of the galactic sig- 
nal in the contaminated region without changing anything in 
the clean region where one wants to perform the CMB analysis. 
The efficiency of the ringing effect removal is shown in Fig.|6l 
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Fig. 6. The error bars ont the power spectrum estimation are 
shown shown in black. The excess power spectrum contribu- 
tion coming from the ringing effect on the Galactic signal is 
shown in blue for a raw Fourier cut on the inital data.. The ef- 
fect of the dust removal by subtracting our estimate of the dust 
signal is shown in green. Adding the smooth sine shape filter 
instead of a raw Fourier cut leads to the red residuals. 

The ringing coming from the Galactic signal is efficiently re- 
duced to a level much smaller than the CMB anisotropics and 
the error bars on the power spectrum estimation. 

4. Conclusions 

We have developped a method that allows an accurate esti- 
mate of the noise power spectrum for CMB experiments ac- 
counting for small but significant effects such as the noise re- 
projection and the pixellisation. This method relies on Monte- 
Carlo simulations and does not require specific platforms to 
be implemented. It is designed to serve as an input to CMB 
power spectrum estimation that generally Requires a precise 
estimate of the noise properties. This method has been suc- 
cesfully applied in the frame of the Archeops data analy- 
sis ( Ben oit et al. 2003al Benoit et al. 2003b) but can be ex- 
tended to any non interferometric experiment. 
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